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We define the "BEC state" to be the many-body wavefunction where all particles are in the same 
one-body state. Using an argument analogous to Anderson's Orthogonality Catastrophe, we argue 
that for interacting particles the amplitude of the BEC state within the many-body wavefunction 
goes to zero in the thermodynamic limit. This does not mean that there is no condensate. However, 
we argue that, if the excitations satisfy the Landau criterion, then the absence of a finite amplitude 
for the BEC state, or the absence of a condensate, do not prevent the definition of the phase function 
(f> from which the superfluid velocity follows. 



I. INTRODUCTION 

It is usually considered that Bose-Einstein conden- 
sation (BEC) is necessary for superfluidity in a one- 
component system^— BEC is invoked because it enables 
one to conceptually label a part of the many-body wave- 
function (the condensate) that can be associated with a 
quantum-mechanical phase (/>, from which the superfluid 
velocity is defined via 

Vs = -W, (1) 

m 

where m is the mass of the particle that constitutes the 
one-component system. Certainly BEC is sufficient for <\> 
to be well-defined, but in the present work we argue that 
BEC does not occur in the thermodynamic limit, when 
the number of particles N — > oo. However, following 
Landau's work,~ which does not comment on this issue, 
we argue that BEC is not necessary for (j> to be well- 
defined. 

In what follows we distinguish between what we call a 
"BEC state" and a "BEC condensate fraction" fo, and 
both are distinct from the superfluid fraction f s & There- 
fore we consider an issue that, to our knowledge, is new 
to the literature. By BEC state we mean a many-body 
wavefunction that consists of a product where all of the 
N particles are in the same one-body state. (Ref. y also 
uses the phrase "BEC state" , which could also mean the 
thermodynamic state of a non-interacting Bose system 
at temperature T — 0.) The conceptual advantage of 
having a finite BEC state is that for this state all of the 
particles have the same phase (up to a constant). There 
can be an infinite number of BEC states, but for a system 
at rest only the zero-momentum state is relevant. 

On the other hand, the "BEC condensate fraction" 
fo is the fraction of particles in the one-body state, av- 
eraged over the many-body wavefunction. Thus, with 
Nq the average occupation number for the many-body 
wavefunction, we have fo — No/N. Hence, in the title, 
by "average BEC" we mean "finite condensate fraction" . 
The one-body state may be "dressed" by an average one- 
body interaction with the other particles; it need not be 
a bare state (such as the zero momentum state). 

This distinction between the BEC state and fo is not 
merely semantic. If each configuration included in the 



many-body wavefunction consists of a product of 80% 
particles in the same one-body state, and 20% particles 
not in this one-body state, then there is a 0.8 conden- 
sate fraction but no BEC because none of the configura- 
tions contains 100% of the particles in the same one-body 
state. In this work we argue that the BEC state ampli- 
tude goes to zero as N -> oo. However, that does not 
imply that fo goes to zero in that limit. In fact, the value 
of fo does not appear in the macroscopic theory of su- 
perfluidity. According to our understanding of Landau's 
viewpoint, the value of fo is irrelevant to our ability to 
define </>. Therefore fo may be zero or it may be non- 
zero, and we can still (as shown in Sec. IIV[) define <fi. In 
fact, for a supersolid, the phase function can depend on 
many-body correlations, although the simplest form of <$> 
is the usual one-body form <fi(r)Z 

Our argument that the amplitude of the BEC state, 
which can be part of a many-body state, has an analog 
in fermion physics. Anderson has given a related argu- 
ment for a non-interacting Fermi system to which a single 
scattering center is added, sometimes called the Orthog- 
onalization Catastrophe^ He shows that, for N fermions, 
each non-interacting one-body wavefunction, when per- 
turbed by the scattering center, has nearly unity overlap 
with the unperturbed scatterer-free states, but the (anti- 
symmetrized) product of the N wavefunctions has nearly 
zero overlap with the product of the N scatterer-free 
states as N — > oo. 9 In the present bosonic case, for the 
unperturbed scatterer-free A^-body state - which is pre- 
cisely the BEC state - the overlap goes to zero even more 
strongly (exponentially with N) than in the fermionic 
case (power law in A^" 1 ). Hence the amplitude of the 
BEC state goes to zero as A" — > oo. 

Sec. |TT] gives a detailed argument why the interactions 
cause the full many-body wavefunction to have essen- 
tially zero amplitude for the BEC state, even when there 
is a finite condensate fraction. Sec. Mil discusses how the 
phase can be defined even when there is no finite am- 
plitude BEC state. Sec. IIVI gives a further discussion of 
the phase. Sec. [V] proposes a test of the idea that there 
is a finite superfluid fraction even when the condensate 
fraction vanishes. Sec. I VII discusses the Landau criterion. 
Sec. IVIll provides a summary and conclusions. 
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II. WHY THE BEC STATE HAS NEGLIGIBLE 
AMPLITUDE 

First consider the non-interacting Bose system at tem- 
perature T = 0. The wavefunction in that case is the 
BEC state. Although it is Bose-condensed, with every 
particle in the same state, as noted by Landau it is un- 
stable to excitations when any object passes through it, 
and therefore it is not a superfluid. Bogoliubov found 
that turning on a weak interaction stabilizes the excita- 
tions, so that the system satisfies the Landau criterion 
for superfluidity. However, interactions also decrease the 
condensate fraction, which if BEC were essential to su- 
perfluidity, would seem to decrease the stability of the su- 
perfluid. Indeed, in principle it should be possible for the 
condensate fraction to be zero (not merely zero amplitude 
of the BEC state), yet if its excitations satisfied the Lan- 
dau criterion the system would be a superfluid. Another 
result of turning on the interaction is that the relative 
parts of the many-body wavefunction, written in terms 
of suitably symmetrized products of non-interacting one- 
body states, have fixed relative phases and fixed ampli- 
tude relative to the BEC state, as discussed in detail in 
Sect. III. 

It might appear that the state found by perturbation 
theory in Bogoliubov's model has a finite BEC state 
amplitude. However, we now argue that is not the 
case. Consider an example in which each of the orig- 
inal one-body states within the many-body state, af- 
ter interactions are included, overlaps by 0.9 with the 
original common non-interacting one-body state. (For 
simplicity it may help to think of the system as having 
only two available one-body states, with amplitudes 0.9 
and y/l - (0.9) 2 ). Then for N particles the many-body 
wavefunction overlaps by (0.9)^ with the original non- 
interacting many-body state, by which we mean the BEC 
state. This overlap vanishes in the thermodynamic limit 
N — > oo. Hence there is no many-body BEC, although 
there is a finite condensate fraction of 0.9. 

This means that, in contrast to the case of a finite am- 
plitude BEC state, there is no obvious label for a state 
whose phase can be used to define a superfluid veloc- 
ity. The viewpoint of Landau, however, is perfectly com- 
patible with the BEC state having negligible amplitude. 
Although Landau does not emphasize the point, the in- 
teractions stabilize the wavefunction in the sense that it 
is a sum of products of one-body terms whose relative 
amplitude and phase are determined by the interactions 
via the Schrodinger equation. (We discuss this in detail 
in Sect. III.) This makes it possible for the system to 
have an overall phase function, which is the requirement 
of Landau's theory of superfluidity. (Landau's theory is 
also subject to the stability of the system relative to a 
moving wall causing spontaneous generation of excita- 
tions.) In principle, many-body effects can modify this 
phase, but they matter only for systems that are not 
translationally invariant^ 

One can argue that it is an oversimplification to begin 



by describing the interacting system description with un- 
perturbed many-body wavefunctions. It would be more 
accurate to consider a Hartree or Hartree-Fock or mean- 
field approximation to the one-body wavefunctions that 
includes an average interaction with the other particles. 
Those wavefunctions would have a condensate fraction 
decreased from unity (e.g., 0.7). If they were exact eigen- 
states the system would still be Bose-condensed with a 
condensate fraction of 0.7. However, going beyond the 
average interaction of mean-field theory, the many-body 
wavefunction becomes more complex on including scat- 
tering processes omitted by the average interaction. If, 
as in the previous example, scattering causes each of the 
mean-field one-body states within the many-body state 
to now overlap, as before, by 0.9 with the common mean- 
field one-body state, then for N particles the many-body 
wavefunction now overlaps by (0.9) with the original 
common non-interacting many-body state. The conden- 
sate fraction would thus vanish in the thermodynamic 
limit. Hence there is no many-body BEC when interac- 
tions are included, although there is a finite condensate 
fraction of 0.7 x 0.9. 

The above argument does not include higher-order cor- 
relations. However, these should be unimportant for 
weak interactions, and we expect they would not invali- 
date the general argument. 



III. GLOBAL PHASE IS DEFINED EVEN 
WITHOUT BEC 

Consider a system of any number of identical Bose 
particles. In the absence of interactions, the eigenstates 
of the system as a whole are products of eigenstates of 
the individual particles. The Landau argument applied 
to the extremely low energy-per-momentum free-particle 
excitations indicates that the system is not superfluid. 5 
Now turn on the interactions. By the fundamental princi- 
ples of quantum mechanics, any eigenstate of this system 
can be written as a sum of configurations, each configura- 
tion being a product state for the non-interacting system. 
By the diagonalization procedure associated with finding 
eigenstates of the Schrodinger equation, each configura- 
tion has a well-defined amplitude and phase relative to 
every other configuration. Only the overall phase of the 
system is unknown. This is true even if the condensate 
fraction fo goes to zero. (It is also true for fermions, and 
for arbitrary mixtures of particles of different types with 
differing statistics.) Thus an overall phase for the system 
can be defined even without a finite fo- The argument 
also applies if there are multiple one-body states with 
macroscopic occupation numbers. 

To be specific, consider conventional perturbation the- 
ory with perturbation V, as applied to the ground state 
'I'j '' of an interacting system I. With 4^ = \ j, nl) de- 
noting the states of the non-interacting system nl, and 
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normalized to unity, we have the well-known result that 

where Ej_„i is the energy of the noninteracting state 
|j, nl) and C is a dimensionless normalization constant. 
This result applies both to a many-body system and 
a one-particle system. Therefore the relative values of 
the amplitude and phase between any two noninteract- 
ing states, mixed in by the interaction, are completely 
determined. 

More generally, even beyond perturbation theory, we 
have 

for some dimensionless C. Hence, even if, as argued 
above, there is no finite BEC amplitude (i.e., the dimen- 
sionless quantity Aq is negligibly small), it is possible to 

attribute an overall phase to , which can be consid- 
ered to be built into C. If the wavefunction is known 
analytically, then one can study whether or not such a 
state is topologically connected] 10 ' 11 If the wavefunction 
is not known analytically, another approach would be to 
employ Quantum Monte Carlo methods to study the su- 
perfluid fraction ! 12 ' 13 For a flowing system the phase can 
be made spatially-dependent, as in the previous section, 
and one can then study whether or not such a state is 
stable according to the Landau criterion. 

How to go beyond a one-body phase, to include higher- 
order correlations, has also been considered. 7 The formal- 
ism for the case of a two-body phase has been developed, 
but not applied, because currently the three-body corre- 
lations needed to compute the two-body phase function 
upper limit are unavailable. Specifically, the wavefunc- 
tion with superflow for i = 1, N particles can take the 
form 

*({!«}) = *o({a:i})exp(t$), (4) 
= ^<l>i(xi) + 2 {N -1) Z^fo' 3 ^) 

i ij 

+ ... (5) 

Here is the ground state wavefunction without flow, 
4>i is the usual one-body phase function familiar from 
liquid 4 He, and (j>2 is a new phase function that includes 
two-body correlations. The above refers to T = 0. 



IV. ON THE PHASE OF THE FIELD 
OPERATOR AND ON ODLRO 

In many cases the one-body phase <f> associated with 
superfluidity is considered to be a consequence of the 
behavior of the two-body correlation function, at a large 



separation. Specifically, it is argued 3, - that one can obtain 
the phase via the one-body density matrix 

{i/y (r')ijj{r)) > noexp[i0(r') — i(fi(r)]. (6) 

| f ' — r | — > oo 

Here ip{r) is the one-body field operator, i/j'(r) is its 
Hcrmitian conjugate, and no is the number density of 
the condensate. This is known as off-diagonal long-range 
orderp^ or ODLRO. We believe that this is an asser- 
tion rather than a proof that no can be finite. Moreover, 
because the wavefunction is defined only up to an over- 
all phase, the above equation is trivial when associated 
with a non-zero no- Because of the Landau argument the 
phase is defined even when no — > 0. 

In some cases one considers the form (ip{f')). Such a 
field-average is in Fock-space, with many different num- 
bers M of particles. If the Fock-space states are uncorre- 
lated, then this field-average must be zeroJ^ If the Fock- 
space states are correlated, then the field-average can be 
non-zero, and is a sum over M, with matrix elements 
having the bra (M + 1| on the left and the ket \M) on 
the right. As a consequence the matrix element clearly 
involves integration over the M variables on the right, 
leaving one variable (call it tm+i) on the left unspeci- 
fied, in addition to the variable r of the field operator. 

We also refer the reader to Ref. 3 for a general discus- 
sion of superfluidity, based on the assumption that there 
is one eigenstate of the one-body density matrix whose 
eigenvalue is macroscopic. This assumption is distinct 
from the assumption that the BEC state has a finite am- 
plitude, but it appears to us that they are likely to be 
equivalent statements. Certainly if there were BEC con- 
densation into a single state, one would expect a single 
eigenstate of the one-body density matrix to have a sin- 
gle macroscopic eigenvalue, and that the corresponding 
eigenfunction would correspond to the BEC state. 

A similar criticism can be made of ODLRO-based argu- 
ments for superfluidity. Nevertheless, because (following 
Landau) the overall phase gradient is well-defined, it is 
unnecessary to use field averages or ODLRO arguments 
to provide a conceptual basis for the phase. 

Finally, we believe there is no compelling reason why 
a system cannot have condensate fraction /o = and yet 
also have excitations that satisfy the Landau criterion. 
In that case one would have a supcrfluid without a con- 
densate. We believe that it is worth investigating this 
possibility in "analytic" many-body theories^ 

V. NUMERICAL TEST FOR BEC AND 
SUPERFLUIDITY 

Ref. H notes that there may be more possibilities for 
a bose fluid than, as seems to be the case for 4 He, so- 
lidification and BEC condensation. Bose systems can be 
simulated numerically, and this suggests a numerical test 
for the necessity of BEC to superfluidity. There exist ac- 
curate analytic theories for liquid 4 He that can be used 



4 



with other potentials, and for which one can expect to 
obtain reasonably accurate results^ If one starts with 
the potential for liquid 4 He and increases its strength we 
expect the condensate fraction to decrease. If one also 
adjusts the shape of the potential, it may be possible to 
find potentials that lead to liquid states (at T = 0) for 
which the condensate fraction is zero. There are then two 
possibilities: either the excitations are stable, or they are 
unstable (in the sense of the Landau criterion). 

Using such a potential in a quantum Monte Carlo al- 
gorithm, the supcrfluid fraction can be computed via the 
winding numberi 12 ' 13 The necessity-of-BEC viewpoint 
predicts that if the condensate fraction is zero, then the 
supcrfluid fraction must be zero. The Landau criterion 
viewpoint predicts that if the excitations are stable then 
the supcrfluid fraction must be non-zero. Such numerical 
studies should be very revealing. 

VI. ON LANDAU'S CRITICAL VELOCITY 
ARGUMENT: CAN WALL EXCITATIONS 
DESTROY SUPERFLOW? 

Ref. [H comments that Landau's critical velocity argu- 
ment most straightforwardly applies only to ions (struc- 
tureless objects) moving through the putative supcrfluid. 
For example, flow past a wall is more complex than that 
of an ion because the wall has excitations, and momen- 
tum and energy transfer from fluid to wall might involve 
the wall's excitations. In practice, vorticity, rather than 
ordinary non-topological excitations, is created at the 
walls, and destroys supernuid flow. 

Nevertheless, consider the possibility that excitations 
of the wall (rather than the fluid) can take up momentum 
and energy, and thereby can destroy superfluidity The 
phonons of the wall material, with a finite velocity, will 
be stable under the Landau criterion. Since phonons are 



the only low-energy excitations for insulators, flow past 
insulators should be stable under the Landau criterion. 

Now consider a confining wall that is a conductor, 
for which quasiparticles are filled to the top of the 
Fermi sea. For flow along z an electron at the Fermi 
level with momentum along x can absorb momentum 
along z, exciting it to above the Fermi level at the cost 
of very little energy. Therefore from energetics alone 
it appears that interactions of He with conducting 
walls can, in principle, destroy superflow. We are 
aware of no evidence for this; perhaps the effect is too 
small to be observed because the electron wavefunctions 
are exponentially suppressed from "leaking" into the 4 He. 



VII. SUMMARY AND CONCLUSION 

We argue that the BEC state amplitude is zero in the 
thermodynamic limit of N — > oo. We also argue that, if 
the excitations of the system satisfy the Landau criterion, 
then a finite one-body condensate fraction is unnecessary 
for superfluidity, and we propose a test of this proposal 
by an analytical-numerical many-body theory (to find a 
liquid state with stable excitations and no condensate) 
followed by a quantum Monte Carlo calculation to deter- 
mine the supernuid fraction. 
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